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We review applications of Dyson-Schwinger equations at nonzero temperature, 
T, and chemical potential, ji, touching topics such as: deconfinement and chiral 
symmetry restoration; the behaviour of bulk thermodynamic quantities; the (T, /i)- 
dependence of hadron properties; and the possibility of diquark condensation. 



1 Introduction 

Confinement and dynamical chiral symmetry breaking (DCSB) are conse- 
quences of the little-understood long-range behaviour of the QCD interaction, 
and developing a better understanding of that behaviour is a primary goal 
of contemporary nuclear physics. It is a prodigious problem whose solution 
admits many complementary strategies. Our approach is to apply a single 
phcnomcnological framework to many observables, thereby identifying the uni- 
fying qualitative features. Non-hadronic electroweak interactions are the best 
observables to study because the probes, the photon and W, Z bosons, are 
very well understood. Following such applications!!! we can infer consequences 
for QCD at extremes of temperature and chemical potential. 

Our tools of choice are the Dyson-Schwinger equations^ (DSEs), which at 
the simplest level provide a means of generating perturbation theory and are 
an invaluable aid in proving renormalisability. However, our interest stems 
from their essentially nonperturbative character. For example, the DSE for 
the quark propagator is the QCD gap equation. Its complete solution contains 
all that is necessary to describe DCSB and yields insights into confinement, 
both of which are absent at any finite order in perturbation theory. Further, 
the Bethe-Salpeter equations (BSEs) are just another form of DSE and these 
equations completely describe meson structure. 

The formulation of the DSEs is straightforward but their solution is not. 
The equation for a particular propagator or vertex (n-point) function involves 
at least one m > n-point function; e.g., the gap equation whose solution is the 
dressed-quark propagator (2-point function) involves the dressed-gluon propa- 
gator, a 2-point function, and the dressed-quark-gluon vertex, a 3-point func- 
tion. Thus in the DSEs we have a countable infinity of coupled equations and 
a tractable problem is only obtained if we truncate the system. This has been 
an impediment to their application: a priori it can be difficult to judge the fi- 
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delity of a particular truncation scheme. However, with expanding community 
involvement this barrier is being overcome as truncation schemes are explored 
and efficacious ones developed. 



2 Gap equation 

The gap equation in QCD is the DSE for the quark propagator: 

Sfip)- 1 ■.= i 1 -pA f {p 2 ) + B f (p 2 )=A f {p 2 )(i 1 -p + M f {p 2 )) (1) 

= Z 2 (i 1 -p + m h f ™) + Z 1 [ g 2 D^(p-q)^- lt2 S f (q)Tl a (q,p),(2) 



■i 



is the dressed-gluon propagator (in Landau gauge, just to be concrete) , (g, p) 



is the dressed-quark-gluon vertex, m 1 ?" 1 is the A-dependent bare /-quark current- 



mass and J A := J A d 4 q/ (2tt) 4: represents mnemonically a translationally- in- 
variant regularisation of the integral, with A the regularisation mass-scale. The 
renormalisation constants for the quark-gluon- vertex, quark wave function and 
mass: Zi(C 2 ,A 2 ), Z 2 (C 2 ,A 2 ) and Z m (C 2 ,A 2 ) := Z 2 (( 2 , A 2 )- 1 Z 4 (C 2 , A 2 ), de- 
pend on the renormalisation point, £, and the regularisation mass-scale. (The 
renormalised current-quark mass is ti/(C) : = ^m lm / m -L 

The qualitative features of the QCD solution of Eq. (||) are known. The chi- 
ral limit is defined by rh = 0, where m is the renormalisation-point-independent 
current-quark mass, and for p 2 > 20 GeV 2 the solution of Eq. (|J) ist 

M Q (p 2 ) ^ - , (4) 



p 2 Qln[pVA 2 QCD ]) 

where j m = 12/(33— 27V/) is the gauge-independent mass anomalous dimension 
and (qq) a is the renormalisation-point-independent vacuum quark condensate. 
The existence of DCSB means that (qq) Q ^ 0, however, its actual value de- 
pends on the long-range behaviour of D^ v (k) and r^ a (q,p), which is modelled 
in contemporary DSE studies. Requiring a good description of light-meson 
observables necessitates (qq}° « —(0.24 GeV) 3 . 

The momentum-dependence in Eq. (Q) is a crucial, model-independent 
result because it is the only behaviour consistent with the definition of the 
vacuum quark condensate as the trace of the chiral- limit quark propagator :H 

,A 

-{qq)° c =N ,lim Z^ 2 , A 2 ) tr D / S (k). (5) 
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Any model that generates 



M (p 2 )~p-' n , »>1 (6) 

will yield (qq)® = from the definition of the quark condensate. 

Confinement is the absence of quark and gluon production thresholds in 
colour-singlet-to-singlet 5-matrix amplitudes. The absence of a Lehmann rep- 
resentation for dressed-quark and -gluon propagators is sufficient to ensure 
thatl3 Therefore the solution of Eq. (Q) can also yield information about con- 
finement, as shown clear lyu for QED3. 

Studies of Eq. (J2j> .lhat employ a dressed-gluon propagator with a strong 
infrared enhancement :lI V(k 2 ) ~ 1/fc 2 , and hence without a Lehmann repre- 
sentation, and Tl a (q,p) regular in the infrared^ yield S(p) that also does not 
have a Lehmann representation. Fine-tuning is not necessary. Such models 
also easily account for D.CSB,0 with_the correct value of (qq)°- 

Contemporary DSEEl and latticel] studies have reopened the possibility that 
V(k 2 ) r~j (k 2 ) p , p < 2, for k 2 ~ 0; i.e., an infrared swpression. The phenomeno- 
logical consequences of this have been re-exploredfl when that V(k 2 ) obtained 
in contemporary lattice simulations is used in Eq. ^) with an infrared-regular 
dressed-quark-gluon vertex, DCSB does not occur and S(p) has a Lehmann 
representation; i.e, there is no signal of confinement. The V(k 2 ) ~ (fc 2 ) p -form 
obtained in DSE studies can be mad£_to support a nonzero condensate via 
Eq. (||), however, its value is typicallyollj only 7-30% of that required to explain 
observed phenomena, and again S(p) does not exhibit signs of confinement .|j 

3 Exploring QCD at nonzero T and fj, 

The dressed-quark propagator at nonzero- (T, fj.) has the general form 

27 • p A(p k ) + i"jiU> k+ C(p k ) + B{p k ) ' 
= -ij ■ pa A (p k ) - i^iUJ k+ o- c {pk) + o- B (Pk) , (8) 

where we have omitted the flavour label, p k — (p,Wk + ), ^fc + = u>k + i[i 
and u>k = (2fc + l)nT, with k £ Z, is the quark's Matsubara frequency. 
The complex scalar functions: A(p,ui k+ ), B(p,u> k+ ) and C(p,u> k+ ) satisfy: 
J-{p,uJk + )* = J~(p, uj-k + -i) i T = A,B,C, and although not explicitly indi- 
cated they are functions only of \p\ 2 and ^ 2 + - The dependence of these func- 
tions on their arguments has important consequences in QCD: it can provide 

a It is difficult to interpret particle-like singularities in coloured Schwinger functions in a 
manner consistent with confinement 
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an understanding of quark confinement and is the reason why bulk thermody- 
namic quantities approach their ultrarelativistic limits slowly. The nonzero- 
(T, fi) gap equation is a straightforward generalisationO of Eq. (g) and the 
Landau gauge dressed-gluon propagator has the general form 

g 2 D^(p, Sl k ) = P^(p, fik) Af® V k ) + A G (p, a k ) , (9) 

{0; /i and/or v = 4, 

%-p,^ = M = l,2,3> (10) 

witn + Pfrv&Pi) = <V - Wi'/tELi PaPa); ii-,v = 1, •••,4, and 

= 2kirT the boson Matsubara frequency. 

In studying the formation of a quark-gluon plasma, two transitions are im- 
portant: deconfinement and chiral symmetry restoration. The simplest order 
parameter for the chiral transition is 

T 

X{t,h) := ReB (p = ) u> ); t := — - 1 , h := — . (11) 

It is a general result that the zeroth Matsubara mode determines the character 
of the chiral phase transition. 

An order parameter for the deconfinement transition is realiseda via the 
Schwinger function: 

00 1 y* 00 

Ab (x,t = 0) := T V dpp sm(px)a Bo {p,u n ), (12) 



n— — 00 



where we have set fi — for illustrative simplicity. If (Jb (p, ui n ) has complex 
conjugate poles, y p , then: 1) it doesn't have a Lehmann representation; and 
2) Ab (i,t = 0) has zeros. The position of the first zero, ^(t), is inversely 
proportional to lm(y p ). Thus 

««>(*):= l/rg^t), (13) 

is a confinement order parameter because «o(i) — > as i — » 0~ indicates that a 
temperature has been reached at which the poles have migrated to the real-p 2 
axis and the propagator has acquired a Lehmann representation. (This order 
parameter can be generalised to qualitatively different functional realisations 
of the absence of a Lehmann representation.) 
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4 Locating the phase boundary in the (T, /i)-plane 

DSE models constrained at T = = fj, can be used to estimate the location 
of the phase boundary. The studies we review all use rainbow truncation: 
Tv(Qui iPui k ) = lv, which is the leading term in a l/iV c -expansion of the vertex; 
and Landau gauge, with a dressed-gluon propagator characterised by 

A F (poJ = V( P n k ;m g ), A (pnJ = V(p Qlt ; 0) , (14) 
V(pn k ;m g ) :=2Tr 2 D^6ok6 3 (p)+V M (pn k ;m g ), (15) 

PQ k = (p,^k), where D is a mass-scale parameter and T>M{pn k ',rag) may be 
large in the vicinity of pf lk =0 but must be finite. 

11 p = 0,T c =? 

The model obtained with D = (8/9)m| and 

16 2 1 - e~ sn fe /(4m ? ) 

2? Mipa k ; m fl ) = — 7T , (16) 

9 s n k 

where sn k := + m 2 g [m 2 , = 8ir 2 T 2 is a gauge boson Debye mass], yields 
a finite-T extension of a nhenomenologically efficacious one-parameter model 
dressed-gluon propagator .□ The mass-scale m* = 0.69 GeV = 1/0.29 fm was 
fixed by requiring a good description of n- and p-meson properties at T = 0. 
At a renormalisation point of £ = 9.47 GeV, m„(£) = 1.1 MeV yields m,r = 
140 MeV. 

This model has coincident chiral symmetry restoration and deconfincmcnt 
transitions at 

T c x = 0.15 GeV = T c Ko (17) 

with mean field critical exponents. Studie^that employ the rainbow truncation 
must give mean field critical exponentsEil because contributions to the gap 
equation that describe the effects of mesonic correlations, which are expected 
to dominate near the transition temperature, can only arise as corrections to 
the vertex. The behaviour of m, and f n is depicted in Fig. 

As a bona fide order parameter f n cx (— t) 1 / 2 , which is illustrated bjz. 
the curve in Fig. |l|. Hence, it follows from the pseudoscalar mass formulae 
f 2 m 2 = 2m u (Q(qq)p that diverges at the critical temperature; i.e., 
TOtt oc (— t)" 1 / 4 , as illustrated. Qualitatively, these two observations indicate 
that at T = T c x there is insufficient attraction in the pseudoscalar channel for 
a bound state to forma and while correlations may persist above these are 
properly identified as a continuum contribution to the pseudoscalar vertex. 
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Figure 1: The pion mass and decay constant are independent of temperature for T < 0.7T*. 

4.2 T = 0,p c =? 

The difficulties encountered in numerical simulations of lattice-QCD at p =/= 
are described in many contributions to this volume. In studies of the gap 
equation it only means that the self energies are complex-valued functions. 
The T = version of the model in the previous section is obtained with 

in Eq. (0). This model hasEl coincident, first order deconfinement and chiral 
symmetry restoring transitions at p c = 0.375 GoV, as measured by the location 
of the zero in the /i-dependent "bag constant" :Ej B(p). It is positive when the 
Nambu-Goldstone phase is dynamically favoured; i.e., has the highest pressure, 
and becomes negative when the Wigner pressure becomes larger, which is why 
p c is the zero of B(p). To gauge the magnitude of p c we note that in a two- 
flavour free-quark gas the baryon number density pb = 2^i 3 /(37r 2 ) so 

H C = 0.375 GeV =*> p B F+dp = 2.9 Po , (19) 

where po = 0.16 fm~ 3 . This may be compared with the central core density of 
a IAMq neutron star: 3.6-4. lpo, while 0.7p c corresponds to po- 

b The calculated value of 6(0) = (0.104 GeV) 4 = 15MeV/fm 3 is similar to that employed in 
bag-like models. 




A-K z mi8\k) 



-[fc 2 /(4m?)] 



(18) 
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In this model decreases slowly as fj, increases, with m 7r (0.7^( c )/m 7r (0) rs 
0.94. At this point begins to increase although, for fi < /i c , m T (/z) does not 
exceed 771^(0), which precludes pion condensation. The behaviour of m w results 
from mutually compensating increases in f% and m(()(qq)^. is insensitive to 
/j until [i w 0.7 fa, when it increases sharply so that /jr(Mc = 0) » 1.25. 

At jtt c , m^r and / w drop discontinuously to zero. The relative insensitivity 
of rri^ and f v to changes in jj,, until very near /z c , mirrors the behaviour of 
these observables at finite-T.H This study reveals an anticorrelation between 
the /i-dependence of / w and that of m f . 



4.3 0, /j^O 

This is a difficult problem and the most complete studies to dateli-M employ 
the simple Ansatz for the dressed-gluon propagator obtained with D 
and T>yi(pn k ', rn g ) = in Eq. (|l5|), and the mass-scale r\ = 1.06 GeV fixedll3 by 
fitting 7T- and p-meson masses at T = 0. With this Unsafe the gap equation is 



S 1 (p,uj k ) = S 1 (p,Cj k ) + -ri 2 j v S(p,U)khu : 



(20) 



and an integral equation is reduced to an algebraic equation whose solution 
exhibits many of the qualitative features of more sophisticated models. 

In the chiral limit Eq. ( pp[ ) reduces to a quadratic equation for B(pk), which 
has two qualitatively distinct solutions. The Nambu-Goldstone solution, with 



B(p k ) 



C{p k ) 




2, 

3" 



otherwise 



Re(pl) < \ 



2r? 2 

1 + -4- , otherwise 

Pi 



(21) 
(22) 



describes a phase of this model in which: 1) chiral symmetry is dynamically 
broken, because one has a nonzero quark mass-function, B(pk), in the absence 
of a current-quark mass; and 2) the dressed-quarks are confined, because the 
propagator described by these functions does not have a Lehmann representa- 
tion. The alternative Wigner solution, for which 



B(pk) = > C(p k )= 1 -(l + Jl+ 2 -f 



(23) 



describes a phase of the model with neither DCSB nor confinement. 
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Here the relative stability of the competing phases is measured by a (T, /in- 
dependent bag constant :L3 B(T,fj,). The line 2?(T, /i) = defines the phase 
boundary, and the deconfinement and chiral symmetry restoration transitions 
are coincident. For \i = the transition is second order and the critical tem- 
perature is T® = 0.159 77 = 0.17 GeV, just 12% larger than the value reported 



in Sec. 4.1. For any // ^ the transition is first-order and the T — critical 



chemical potential is /i° = 0.3 GeV, « 30% smaller than the result in Sec. 4.2 
The_quark pressure, P q , is easily calculated and P q = in the confined 
domain.Ej However, this does not mean that the vacuum is unaffected by 
changes in (T,fi). On the contrary; e.g., in the models described above, the 
condensate evolves with these changes, as it must because it is a dynamical 
quantity. At each (T, jj) the properties of the hadronic excitations are calcu- 
lated in the evolved vacuum and the modification of the quark-constituents' 
propagation characteristics, which the condensate's modification represents, 
makes a significant contribution to the (T, /^-dependence of those properties. 

In the deconfined domain, P q slowly approaches the ultrarelativistic, free 
particle limit, Pur, at large values of (T, /t); e.g., at T ~ 0.3 77 ~ 2T°, or 
^ r*j -q ~ 3/4;, P q ~ 0.5 Pur- This behaviour results from the persistence of mo- 
mentum dependent modifications of the quark propagator into the deconfined 
domain, as evidenced by C ^ 1 in Eq. fl23]), which also entails a "mirroring" of 
finite-T behaviour in the /i-dependence of the bulk thermodynamic quantities. 

The (T, /^-dependence of vacuum and meson properties is easily calculated 
in this model; e.g., the vacuum quark condensate is 

In 



(W)=^ 3 ^E / ' dyy 2 Re(J\-yZ-u?) (24) 



1=0 



T = T/rj, [i = /i/77; Imax is the largest value of I for which uif < (1/4) + /i 2 
and this also specifies w/ max , A 2 = Ojf — ujf, pi = (y, uji + ip). At T = = fj,, 
(-(qq)) = ?7 3 /(807r 2 ) = (0.11 ry) 3 . "Obvious from Eq. (||) is that (-(qq)) 
decreases continuously to zero with T but increases wilh_ji, up to ^ C {T) when 
it drops discontinuously to zero: as observed elsewhere. clila That behaviour is a 
necessary consequence of the momentum-dependence of the quark self energy, 
with the finite- (T, fi) behaviour of observables determined by 

Re(o; 2 u] ) d ~[7r 2 T 2 -M 2 ] d , (25) 

where d is the observable's mass-dimension. This is confirmed in the chiral 
limit expression 



In 



fl = V 2 ^f g f (l + 4 f - 4^ 2 - § A 2 ) . (26) 
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The anticipated combination fi 2 — ujf appears and even without calculation it 
is clear that / ff will decrease with T and increase with /i. 

The (T, p)-response of meson masses is determined by the ladder BSE 

T M {Pk\ Pi) = ~ Re { 7m + iPi) r M (p i; Pg) - ~P*) 7 M } , (27) 

where P^ := (P, f^), with the bound state mass obtained by considering P^=o. 
In this truncation the lo- and p-mesons are degenerate. 

The pion solution of this equation is T 7T (P ) — 75(16*1 + 7 • P 9 2 ) and, con- 
sistent with what we saw_a.bove, the mass is (T, ^-independent, until very near 
the transition boundaryli3 For the p-meson the solution has two components: 
one longitudinal, 9 p+ , and one transverse, Q p _, to P. Equation (27) yields an 
eigenvalue equation for the bound state mass, M p ±, and using the chiral-limit 
solutions, Eq. (|2l|), one finds immediately that 

M 2 _ = r] 2 /2, independent of T and //. (28) 

Even for m ^ 0, M p - changes by < 1% as (T, p) are increased from zero 
toward their critical values. This insensitivity is consistent with the absence of 
a constant mass-shift in the transverse polarisation tensor for a gauge-boson. 
For the longitudinal component one obtains in the chiral limit: 

f2 _ 1 2 a(.,2 „2r^ 
2 

The combination /1 2 — ir 2 T 2 again indicates the anticorrelation between the 
response of M p+ to T and its response to p, and, like a gauge-boson Debye 
mass, that M p+ rises linearly with T 2 for fi = 0. The m ^ solution for the 
longitudinal component is semiquantitatively the same. 

The BSE yields qualitatively the same behaviour for the 0-meson. The 
transverse component is insensitive to T and fi, and the longitudinal mass, 
M0_|_, increases with T and decreases with //. Using r\ = 1.06 GeV, Af^± = 
1.02 GeV for m s = 180 MeV at T = = p. 

In a 2-flavour, free-quark gas at T = nuclear matter density corresponds 
to fjb = Ho := 260 MeV = 0.245 77 and the algebraic model yields 

Mp+ifio) ~ 0.75M p+ (n = 0) , M^+ifio) ~ 0.85M^+(p = 0) . (30) 



A/ p 2 + = V - 4 (p 2 - n 2 T 2 ) . (29) 



Section 4^2 indicates that a better representation of the ultraviolet behaviour 
of D pu (k) increases the critical chemical potential by 25%. This suggests that 
a more realistic estimate is obtained by evaluating the mass at fjf = 0.20 77, 
which yields 

M p+ (n' ) « 0.85Af p+ (p = 0) , A^+^o) « 0.90M*+(/i = 0) ; (31) 
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a small, quantitative modification. The difference between Eqs. (|30|) and (31) is 



a measure of the theoretical uncertainty in the estimates in each case. Pursuing 
this suggestion further, [i = y/2 (i' Q , corresponds to 2po, at which point M u+ = 
M p+ w 0.72 M p+ (fi = 0) and M 0+ « 0MM^+(jJt = 0), while at the T = 



critical chemical potential, which corresponds to approximately 3po in Sec. 4.2 
M u+ = M p+ w 0.65 M p+ (^ = 0) and M 0+ w 0.80 M^+(/i = 0). These are the 
maximum possible reductions in the meson masses. 

5 Diquark condensation 

A direct means of exploring the possibility that SU(iV c ) gauge theories might 
support scalar diquark condensation is to study the gap equation satisfied by 

-{-a>(p)\\t] 15 c(S(-p)-T&)> (32) 

where Sip)- 1 = i<y ■ p A{p 2 ) + B{p 2 ), {\\ 7 i = l...d c ,d c = N C (N C - l)/2} are 
the antisymmetric generators of SU(JV C ), C — 7274 is the charge conjugation 
matrix, and here we consider SU(iVy = 2). A' l (p) ^ for any i indicates the 
formation of a diquark condensate. Sip)- 1 is a matrix in the space of quark 
bispinors: 

Q & = 7i (S) ' Q(x) = {q(x) g_c(x)) ' (33) 

q c := —qC, q c :— q T C. The gap equation is 

W-W+^cgfao), (34) 

where <5>o(p) _1 = diag(i7 ■ p + m, C(— ■ p + m) T C^) and the form of £(p)jj 
specifies the theory and its truncation. This approach avoids a truncated 
bosonisation, which in all but the simplest models is a procedure difficult to 
improve systematically and prone to yielding misleading results. 

SU(-/V C = 2): Using the rainbow truncation and a Feynman-like gauge for illus- 
trative simplicity, the T = = fx gap equation in this theory yields 

pW)-l) = f [ \ 2 D{p-k)p-k^±, (35) 
B{p 2 ) -m = \z\J 9 2 D(p k) (36) 
^ 2 )=\E\j\ 2 D{p-k)^, (37) 
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where d(p 2 ) = p 2 A(p 2 ) 2 + B(p 2 ) 2 + A(p 2 ) 2 , and the pseudo-reality of SU(2) is 
responsible for the identical couplings in Eqs. ( p6| ) and (|37|). Clearly, for m = 
the theory admits degenerate and indistinguishable quark and (colour-singlet) 
diquark condensates. This result is valid independent of the truncation and 
gauge, and is just one of the manifestations of Pauli-Giirsey symmetry. SimL. 
larly, mesons and baryons, which are diquarks in SU(iV c = 2), are degenerate.E3 
The phase structure of this theory at nonzero-(T, //) can certainly be very rich. 

SU(iV c = 3): In this case the interesting possibility is the existence of a colour 
antitriplet diquark condensate: {A A }i=i,2,3 = {A 2 , A 5 , A 7 }. Choosing the con- 
densate to point in the A^-direction, the bispinor propagator separates into two 
pieces, one parallel and the other perpendicular to the condensate's direction: 



s(pY 



Here, since 





f 


Snip)- 1 it o 







1 


o s^p)- 1 












-A*(p)A^76 


S||(P)"^I 












S A 


-(f)- 1 



(38) 



ST 



= h 



3 II 



V3 



S7 1 - 



(39) 



the A a 5A a interaction in the gap equation couples the parallel and perpendic- 
ular components. In rainbow truncation and using a Feynman-like gauge, the 
gap equation yields 



P 2 {A^p 2 )-l) 
p 2 (A ± (p 2 )-l) 
B\\ (p 2 ) — m 
B ± (p 2 )-m 
AV) 



g 2 D(p — k) p ■ k 



A ± (k 2 ) , 5 ^||(fc 2 ) 



9 2 D( P 



9 2 D(p 



k)p ■ k 



dAk 2 ) 
2A ± (k 2 ) 



3 d\\(k 2 
, Mk 2 ) 



k) 

g 2 D{p-k) 
\ J A g 2 D(p-k) 



3d ± (k 2 ) d\\{k 2 ) 

n B ± (k 2 ) i i i?ii(fc 2 y 

d^k 2 ) 
4B ± (k 2 ) 



3 d\\(k 2 ) 



3 d±(k 2 ) 
A\k 2 ) 



d\\{k 2 ) 



d\\{k 2 ) ' 

d ll (p 2 )=p 2 A ll (p 2 ) 2 + B ll (p 2 ) 2 + (A 1 (p 2 )) 2 and d ± (p 2 ) = p 2 A ± {p 2 
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(40) 

. (41) 
(42) 

(43) 
(44) 
B ± (p 2 ) 2 . 



The class of models hitherto applied in exploring diquark condensational 
can be characterised as those in which J (Kl^p ■ k D(p — k) = 0. In this class 
A|| = A± = 1 and when B\\ = B± = 0, which is always a solution, A 1 ^ if 
the coupling is large enough. Hence such restricted models admit a rich phase 
structure at nonzero-(T, /i) because the || <->_L coupling is eliminated. 

However, if one includes the next order contribution to the kernel of the 
gap equation, the picture can change. One studytj suggests that in that case, 
even without the |) ^->_L coupling, A 1 = is the only solution at T = = /J,. 
The effect at nonzero- (T, /z) of correcting the kernel has yet to be investigated 
but this result signals a need for caution in making inferences about the phase 
structure of QCD based on the rainbow-like truncation of this class of models. 

The more general class of models in which J dVtjp ■ k D(p — k) ^ can be 



exemplified by the confining model introduced in Sec. 4.3. In that case, if we 
consider Bu = B± = 0, the gap equation is solved with 

P %( P 2 ) 2 + (AV)) 2 = ^ 2 (45) 

and (setting 77 2 — * 1) 

A \\ (P 2 ) = 1(7 + 3^9 + W) > A±( P 2 ) = \ (3 + ^9 + 2/p 2 ) . (46) 

However, inserting the result for A\\(p 2 ) into Eq. ( |4^ ) yields (A 1 (p 2 )) 2 < 
for all p 2 ] i.e., A 1 ^ 2 ) = 0, even in the rainbow truncation. Thus diquark 
condensation at T = = /1 is blocked by the |j couplingj—. We expect 
that this conclusion will be reinforced if the kernel is improved.! 2 !! The effect 
of /i 7^ has not yet been explored but this result too advocates caution 
in making inferences about the phase structure of QCD based on the simple 
models hitherto employed. 



6 Epilogue 

Hadron observables are insensitive to the behaviour of the interaction at p 2 < 
Aq CD and the rainbow truncation of the gap equation is quantitatively reliable 
for p 2 > 1 GeV 2 . Thus the model dependence in our approach is contained in 
an apparently small domain. However, as illustrated by the material presented 
in this volume, even that small domain of uncertainty admits a large variety 
of possibilities; although apparently distinct Ansatze may really be different 
realisations of the same phenomena. It is also a crucial domain, covering 
that in which mesonic correlations (vertex corrections) can influence quark 
propagation characteristics, an effect that may become qualitatively important 
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as the phase boundary is approached. Much has been achieved in localising the 
model-dependence but more must be done to further ameliorate it. Discussions 
of the type represented by this volume are crucial to that endeavour. 
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